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Abstract 

We study global dynamics of the Poisson-Nernst-Planck (PNP) system for flows of two types 
of ions through a narrow tubular-like membrane channel. As the radius of the cross-section of the 
three-dimensional tubular-like membrane channel approaches zero, a one-dimensional limiting 
PNP system is derived. This one-dimensional limiting system differs from previous studied 
one-dimensional PNP systems in that it encodes the defining geometry of the three-dimensional 
membrane channel. To justify this limiting process, we show that the global attractors of the 
three-dimensional PNP systems are upper semi-continuous to that of the limiting PNP system. 
We then examine the dynamics of the one-dimensional limiting PNP system. For large Debye 
number, the steady-state of the one-dimensional limiting PNP system is completed analyzed 
using the geometric singular perturbation theory. For a special case, an entropy-type Lyapunov 
functional is constructed to show the global, asymptotic stability of the steady-state. 
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1 Introduction 



Poisson-Nernst-Planck (PNP) systems serve as basic electro-diffusion equations modeling, for ex- 
ample, ion flow through membrane channels, transport of holes and electrons in semiconductors 
(see, e.g., [Tj l2j |2T| I27j and the references therein). There are many excellent works on derivation 
of PNP systems from Boltzmann equations (see PQ and reference therein) by assuming the colli- 
sion time is much smaller than the characteristic time. The PNP systems have been studied under 
various physical relevant boundary conditions such as the non-flux, homogeneous Dirichlet and Neu- 
mann boundary conditions. For those types of boundary conditions, in addition to the total charge 
conservation and the existence of various first integrals, Boltzmann H-functionals or entropy- like 
functionals are successfully constructed, which, together with the advances of Csiszar-Kullback- 
type or logarithmic Soblev inequalities, are applied to investigate the asymptotic behavior of the 
PNP systems and stability of steady-state or self-similar solutions (see, e.g., [3lHl[8ll9l [T0l[Ml[25] ). 
In the context of ion flow through membrane channels, it is physically unreasonable to impose the 
above mentioned boundary conditions on the whole boundary, particularly at the two "ends" of the 
channels. Instead, non-homogeneous Dirichlet conditions on the two "ends" are typically assumed. 
PNP systems supplemented with this type of boundary conditions result in quite different dynam- 
ical behavior. The total charges within the channels are not conserved and entropy-like functionals 
are not available in general, and, most significantly, the asymptotic behavior is different as showed 
in special cases in the present paper. 

In this work, we start a systematic study of the PNP systems modeling ion flows through narrow 
tubular-like membrane channels in physiology. For definiteness, consider flow of two types of ions, 
Si and S2, one with the positive valence a\ > and the other with the negative valence — a% < 0, 
passing through a membrane channel with length normalized from X = to X = 1. Denote the 
concentrations of S\ and S2 at location (X,Y,Z) and time t by c\(t,X,Y, Z) and C2(t,X,Y,Z). 
Then the electric potential $(t, X, Y, Z) in the channel is governed by the Poisson equation 

A$ = -A(«ici - a 2 c 2 ), 

where the parameter A is the Debye number related to the ratio of the Debye length to a character- 
istic length scale. The flux densities, J\ and J2, of the two ions contributed from the concentration 
gradients of the two ions and the electric field satisfy the Nernst-Planck equations 

Di(Vci + aiCiV$) = -Ji and D 2 (Vc 2 - a 2 c 2 V$) = -J 2 , 

and the continuity equations 



where D\ and Z? 2 are the diffusion constants of ions Si and S2 relative to the membrane channel. 
The Poisson-Nernst-Planck system is thus given by 

A3> = - \(a\ci - a 2 c 2 ), 

_L =D X V- (Vci + aiciV$), (L1) 

—± =D 2 V ■ (Vc 2 - a 2 c 2 V$). 
at 

PNP systems have been studied by many authors (see, e.g., [Tl[2l[3lBl[8l[9l[T0l[T5l[T6l[T71 
EU [231 [25l [271 [28]). Many works have been attributed to the one-dimensional PNP systems and 
particularly the steady-state problems (see, e.g., pj)J [HI [23 EZ] ) • Consideration of one-dimensional 
PNP systems is motivated naturally by the fact that membrane channels are narrow. To make this 
reduction more rigorous, we present a mathematical framework based on the ideas in |12l [13] and 
investigate mathematically the limiting process as the three-dimensional domain shrinks to a line 
segment. More specifically, starting with the situation that J7 e is a revolution domain about its 
length (e is related to the maximal radius of the cross-sections of the channel) , we will derive a one- 
dimensional limiting PNP system as e — > 0. Differing from the simple one-dimensional version of 
the PNP system, this limiting PNP system encodes the defining geometry of the three-dimensional 
channel. As the first step in justifying the limiting process, we show the upper semi-continuity of the 
global attr actors A 6 of the three-dimensional systems at e = 0. The existence of global attractors 
for the PNP systems can be found in [9j. It is expected that if the one-dimensional limiting 
system is structurally stable, then its dynamics determines that of three-dimensional system for 
e > small. We will thus examine, in this paper, the steady-state problem of the one-dimensional 
limiting system. For large Debye number, the steady-state problem can be viewed as a singularly 
perturbed one. We show that this problem can be completely analyzed using the geometric singular 
perturbation theory as in [23] . 

The rest of the paper is organized as follows. In Section 2, we give detailed formulation of our 
problem. The domain for the three-dimensional PNP system will be specified and, as the domain 
shrinks to a one-dimensional segment, a one-dimensional limiting PNP system is derived. We then 
state our results on the upper semi-continuity of attractors, on the singular boundary value problem 
of the steady-state PNP system. The proofs are provided in Sections 3 and 4. In Section 3, after 
some technical preparations, we show the upper semi-continuity of attractors. Section 4 is devoted 
to the geometric analysis of the singularly perturbed steady-state problem of the one-dimensional 
limiting PNP system. At the end of this section, a special case is studied for which a if-function 
is found and used to establish the asymptotic stability of the steady-state. 
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2 Formulation of the problem and the statements of main results 
2.1 Three-dimensional PNP and a one-dimensional limit 



We start with setting up our problem. The membrane channel considered here is special and will 
be viewed as a tubular-like domain fi e in M 3 as follows: 

n e = {(x, y, z) o < x < l, y 2 + z 2 < g 2 (x, e )}, 

where g is a smooth (at least C 3 ) function satisfying 

g(X,0) = and g (X) = ^-(X,0)>0 for X G [0,1]. (2.1) 

ae 

The positive parameter e measures the sizes of cross-sections of the membrane channel. For a 
technical reason (used in Lemma l3.ip . we also assume that 

*(M-*<i,«)-o. 

The boundary <9f2 e of O e will be divided into three portions as follows: 



L t ={(X,Y,Z) G dn e 

k e ={(x,y,z) g an £ 
m £ ={(x,y,z) g an £ 



x = o}, 
x = iy, 

Y 2 + Z 2 = g 2 (X,e)}. 



Thus, L t and R e are viewed as the two ends of the channel and M e the wall of the channel. The 
boundary conditions considered in this paper are 



$Il 6 = 0o>O, $|£ £ =0, c k \ Le =l k >0, Ck \ Ae =r k >0, 

d£ _ Oct _ ( 2 - 2 ) 



where 0q, l k and (A; = 1,2) are constants, and n is the outward unit normal vector to M e 
Although the most natural boundary conditions on M e would be the non-flux one 

dcx d§\ , (dc 2 d§\ , 

~7\ Haici— L> = — a 2 c 2 ^— U> =0, 



an an / Me \on on J Me 

the above homogeneous Neumann conditions on M e are reasonable (they are the consequences of 
the non-flux and zero-outward electric field conditions). 

In this paper, we are interested in the limiting behavior of the PNP system when the three- 
dimensional tubular-like domain Q e collapses to a one-dimensional interval as e — > 0. Naturally we 
expect a one-dimensional limiting system whose global dynamics is comparable with those of PNP 
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systems for e > small. This important idea was applied by many researchers in studying the 
dynamics of equations defined on thin domains (see, e.g., [12] CC2 [26], [29] ) . We follow the procedure 
in [13] to derive a one-dimensional limiting system but avoid expressing differential operators and 
transformations in local coordinates. As a result, the expected one-dimensional limiting system is 
more transparent. 

To derive the limiting PNP system, we transfer the e-dependent domain f2 e into a fixed domain 
0, = [0, 1] x D, where B is the unit disk, by applying the following change of coordinates: 

Y Z 



x = X, y 



g(X,e)' g(X,eY 



(2.3) 



In the sequel, we denote by L, R and M, respectively, the boundaries of £1 corresponding to 
L e , R e and M e under the transformation. Let J denote the Jacobian matrix of the change of 
coordinates. Then, 



J 



d(x, y, z) 



1 



d(X,Y,Z) g* 
with det(J _1 ) = g 2 (x,e), and 



g 2 o o 
-ggxy g o 
-gg x z o g 



r 



- g A 



g 4 -g 3 g x y 
-g 3 gxy g 2 + g 2 g 2 x y 2 
-g 3 g x z g 2 g 2 x yz 



d(X, Y, Z) 

d(x,y,z) 



-g 3 g x z 
g 2 g 2 x yz 
g 2 + g 2 g 2 z 2 



l o o 

g x y g o 
g x z o g 



The following result, which can be verified by direct computations, is useful for a clean derivation 
of a limiting PNP system. 

— ► W 1 , ip(p) = q, be a diffeomorphism, and let J(q) = f 2 (V , ~ 1 (?)) be the 



Lemma 2.1. Let ip : R n - 

Jacobian matrix and d(q) = (det J(q))^ 1 ■ If a(p) = (3(ip{p)) : 
the gradients in the two coordinates are related as 



is a smooth function, then 



Further, if £ £ (%)§g 
i=i 

then F(p) = f(tp(p)) satisfies 



V p a(p) = J T {q)V q P{q)- 
for all i = 1, • • • , n, and f : 



is a smooth vector field, 



V P • F(p) 



1 



d{q) • " 

and hence, the Laplace operators are related as 

1 



V,- (%)J(<z)/(g)) 



A p a(p) 



d(q) 



V q ■ (d{q)J{q)J T {q)V q f3{q)) ■ 
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It can be checked that the change of variables in (|2.3p with p = (X, Y, Z) and q = (x,y,z) satisfies 

Therefore, applying Lemma l2.lt system (jl.ip can be rewritten, in terms of (x,y,z), as follows. 

V • (g 2 JJ T V$) = -A(aici - a 2 c 2 ), 

9 

^ = 5r V • (5 2 ^ T Vci + aici 5 2 JJ r V$), (2.4) 



dt g 
dc 2 _ D 2 
~dt~~g 



2 V • (g 2 JJ T Vc 2 - a 2 c 2 g 2 JJ T V$), 



with the boundary conditions 

®\l = 4*o, $\r = 0, c k \ L = l k , c k \ R = r k , 

2.5 

(W, JJV)| M = (Vc fe , JJV)| M = 0, 

where fc = 1, 2 and v is the outward unit normal vector to M. 

By inspecting the structural dependence of JJ T on e, we expect the one-dimensional limiting 
PNP system to be 



7 ^{ 92 °^) = ~ x{aiCl ~ a2C2) 

dc\ D\ d ( 2 9 2 d 



dc 2 _ D 2 d f 2 d 2 d 

<9i <7q cte \ Q dx° 2 a2C2 ^° dx 

on x S (0, 1) with the boundary conditions 

$(i,O) = o , *(M) = 0, c fc (t,0)=/ fc , c fc (t,l)=r fc , (2.7) 

where go(x) is defined in fl2.1|> . 

It was shown in [9] that, for any e > 0, the three-dimensional PNP system has a global attractor 
A e which is a compact subset and attracts all solutions with respect to the norm topology of H 1 x H 1 . 
This result is based on an invariant principle discovered in |10l [TTT [9J [28] for the van Roosbroeck 
models of semi-conductor. The PNP systems are basically the same as the van Roosbroeck models 
and we recall the invariant principle using the above setting. 

Proposition 2.2. Let M be a positive constant with 

M > maxjaiii, a%ri, a 2 l 2 , a 2 r 2 }, 
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and let S be the subset o/iJ (O e ) x ^(Q^ given by 

£ = {(ci,c 2 ) G F 1 ^) x F 1 ^) : < ma < M, < a 2 c 2 < M}. 

Then £ is positively invariant for the PNP system. More precisely, if the initial datum (ci(0), 02(0)) 6 
£ and (ci,C2) is i/ie solution of the PNP system, then (ci(t),c 2 (t)) E £ /or a// i > 0. 

We remark that, for PNP systems with three or more types of ions, the above invariant principle 
is not available. It is not clear to us whether or not a similar principle still holds in this case. PNP 
systems with more than two types of ions are worth further studying. 

The results in [1U\ \TT\ [22] show also that the one-dimensional problem (|2.6|) - (|2.7p has a posi- 
tively invariant region 

S = {(ci,c 2 ) € ^(0,1) x H\0, 1) : < aid < M, < a 2 c 2 < M}, 

where M is the constant in Proposition 12.21 an d problem (|2,6p - (|2,7p is globally well-posed in So 
and has a global attractor Aq in £ H F 1 (0, 1) x F 1 (0, 1). 

Our first result claims that the global attractors A e of the three-dimensional PNP systems are 
upper semi-continuous to the global attractor Aq of the one-dimensional limiting system as e —* 0, 
which partially justify the limiting process. 

Theorem 2.3. The global attractors A e of the three-dimensional PNP systems are upper semi- 
continuous at e = 0, that is, for any r\ > 0, there exists a positive number ei = ei(r/) such that for 
all < e < ei and all w 6 A e , 

distx £ (w,A ) <n, 

where X e = {w : \\w\\ 2 Xc = |M|| 2 + ||Viu||| 2 + ^rHu^H^ + ^11^11^2 < °o}. 

2.2 Steady-state problem of the one-dimensional limiting PNP system 

The steady-state of problem ()2.6p and (|2.7p can be rewritten as 



(2.8 



V 2 -T ( K x )^r- ) = -h{x){aci - (3c 2 ), ^ = 0, ^- = 0, 

dx \ dx J dx dx 

7 i \ dc\ , , \d(j> T , , ,dc 2 , .dd> 

~dx + OCl dx = ~ ~dx ~ ~d~x = ~ 

where fj, 2 = 1/A, J\ = J\/D\, J 2 = J 2 /D 2 and h(x) = g^x), and the boundary conditions are 

0(0) =4>o, c 1 (0) = l 1 , c 2 (0)=h, 
0(1) = 0, ci(l) = n, c 2 (l)=r 2 . 
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Since A is large, we can treat the problem (|2.8p and (|2.9p as a singularly perturbed problem with 
/i as the singular parameter. We will recast the singularly perturbed PNP system into a system of 
first order equations. 

Denote derivatives with respect to x by overdot and introduce 

r = x, u = fih(T)<fi, v = — h(r)(aicx — CK2C2), and w = a\ c\ + a 2 c 2 . 
System f|2.8j) becomes 

1 1 ■ Mt) , T 

[10 =-r--—u, au = v, av = uw + a — — v + n(axJ\ ~ «2 
/i(r) /i(r) 

. ai«2 «2 - «i M , 2 t , 2 T \ (2.10) 

n/(r) n(r) h(r) 

j\ =0, j 2 = 0, t = 1. 
System (|2.10p - i/ie sZow system - will be treated as a dynamical system with the phase space 

M 7 = {(<j>,U,V,W,Jl,J 2 ,T)} 

and the independent variable x will be viewed as time. The boundary condition (|2.9p becomes 

=0o, u(0) = -/i(0)(ai/i - a 2 l 2 ), w(0) = a 2 L x + (3 2 L 2 , r(0) = 0, 
cj){l) =0, u(l) = -ft(l)(ain - a 2 r 2 ), = afn + a 2 r 2 , r(l) = 1. 

Setting (j, = in system (|2.10p . we get the limiting slow system 

1 

= u, = v, = uw, 
h{T) 

aia 2 a 2 -a>! ( 2 . 12) 

ii =0, j 2 = 0, f = 1. 

The set = = v = 0} is called the slow manifold which supports the regular layer of the 
boundary value problem. The regular layer will not satisfy all conditions in (|2,lip if a 2 l 2 — a\l\ 7^ 
or a 2 r 2 — a\r% 7^ 0, and this defect has to be remedied by boundary layers. The boundary layer 
behavior will be determined by the fast system resulting from the slow system (I2.10p by the rescaling 
of time x = Thus, in terms of £, the fast system of (|2.10p is 

' — u, u = v, v = uw + ji— — — -v + /u(«i J\ — a 2 J 2 ), 



(2.11) 



h(r) ' ' " h(r 

a\a 2 a 2 — a.\ 

J[ =0, J' 2 = 0, t' = h 



I «1«2 «2 ~ «1 M / 2 T , 2 t \ (2.13) 

w = u?r \ uv ~ { r? \ '"•w - — — -(a! Ji + a 2 J 2 ), v 7 

rr(r) n(r) n(r) 
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where prime denotes the derivative with respect to the variable £. The limiting fast system at fi = 
is 

(/ 1 / / / aia 2 a 2 — ai 

-U, U = V, V = UW, W = T7T7 — rUV H rT"^ — W, 



J[ =0, 4 = 0, r' = 0. 

The slow manifold Zq is precisely the set of equilibria of (|2,14p . 

Concerning the steady-state problem of the one-dimensional limiting PNP system, we have 



Theorem 2.4. Assume that a\l\ ^ a 2 h and a\r\ ^ a 2 r2 (otherwise, see Remark \4-1\ ). For \i > 
small, the boundary value problem 112. 10\) and h2.11\) has a unique solution near a singular orbit. 
The singular orbit is the union of two fast orbits of system \2.1J$ representing the boundary layers 
and one slow orbit of a blow-up system of \2.12\) (see Section 4 for details) for the regular layer. 
The limiting flux densities are explicitly given by 

( \ ( "2 "1 "2 "1 \ 

Di (In j 1 — ai<po) i {aih) 011+012 {a2h) ai+a2 - {am)^ 2 (a 2 r 2 ) a ^) 

J\ =D\Ji = 7 2 \ — ; > 

I «it»2 i n J ^i—lnQ] f 1 h- 1 (x)dx 

/ \ / a 2 "1 "2 °1 

D2 ( In ^ + a2(po ) ( (ai/i) Q1+a2 (a 2 l 2 ) ai+a2 — ("in) " 1+a2 (a 2 r 2 ) 011+0,2 

J 2 =DoJ: 



2-^2 = T-jT 

2 



a 1 +a 2 I 



Remark 2.1. Note that the factor J Q h~ 1 (x)dx on the denominators in the expressions for the 
flux densities J\ and J2 reflects the effect of the geometry of the three-dimensional channel. Let's 
compare this effect with that of a cylindrical channel where the wall is defined by {Y 2 + Z 2 = e}. 
In this case, the corresponding integral factor on the denominators for the flux densities J\ and 
J 2 is 1. The volume of the channel is ire 2 . For general channels that we considered here, we thus 
assume 

Vol = / irg 2 (x; e) dx = ire 2 . 
Jo 

From which we have Jq h(x) dx = 1. Therefore, 

l=( [ hr 1 / 2 {x)h 1 / 2 {x)dx] < [ h- 1 (x)dx [ h(x)dx = [ h~ 1 (x)dx. 
\Jo J Jo Jo Jo 

The inequality indicates that the more complicated the geometry of the channel, the smaller the 
flux for the ion flow, which agrees with our common sense. 
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3 Upper semi-continuity of attractors 



3.1 Homogenization of boundary conditions 

In this section, we convert the non-homogeneous Dirichlet boundary conditions on L U R in (|2,5p 
to homogeneous ones, while keeping the homogeneous Neumann boundary conditions on M. For 
this purpose, the following technical result is needed. 

Lemma 3.1. Let h : [0, 1] — > R be a smooth function. Then, for any e > 0, there is a func- 
tion H e : n e -» K such that # e pf,0,0) = h{X), H e {<d,Y,Z) = h{0), H e (l,Y,Z) = h(l), and 
(VH"(X, Y, Z), n) = for (X, Y, Z) G M e . 



Proof. We provide a specific construction of a function For convenience, hereafter, we denote 



by flf'(jsr, e) = M(X, e). For any e > and X G [0, 1], let X = ij) e (t, X ) be the solution of 



dX _ J(X,e) 

dt ~ g(X,e) [ } 

with ip e (0, Xq) = Xq. It is easy to see that ip e (t, Xq) is even in t from the equation. Since g'(0, e) = 
g'(l,e) = 0, i) e (t,0) =0and^ e (t,l) = 1 for all t. Therefore, for any (X, t) G [0, l]x[0,g(X,e)], there 
is a unique Xq G [0, 1] such that X = ^> e (i, -Xo), and hence, for any (X, Y, Z) G £l t , there is a unique 
X Q G [0, 1] such that X = ^ € (VY 2 + Z 2 , X ). Set H e (X, Y, Z) = h(X ) if X = ip e (\/Y 2 + Z 2 , X ). 
Then, iT e (X,0,0) = h(X), H e (0,Y,Z) = h(0) and H e (l,Y,Z) = h(l). It remains to show that, for 
(X, Y, Z) G M e , (VH e (X, Y, Z), n) = 0. For any X G [0, 1], the set 



D{Xq) = {(X,Y,Z) : X = ^(Vy 2 + Z 2 ,X )} = {(X,Y,Z) : H(X,Y,Z) = h(X )}, 

is a level set of H e . Note also that the curve {(X, Y,0) : X = ip e (Y, Xq)} lies on D{Xq) and it is a 
solution curve to (|3.ip if Y is viewed as the t- variable. Therefore, at (X,Y,0) = (X,g(X,e),0) G 
D(X )nM € , the vector 

-Y 9 -^-,l,o) =(-g'(X,e),l,0) 

is tangent to D(Xq), and hence, (S/H e (X, g(X, e), 0), (— g'(X, e), 1, 0)) = 0. Since n is parallel to 
(-g'(X,e),l,Q), {VH e (X,g(X,e),0),n) = 0. Due to the rotation symmetry of M e and H € about 
the X-axis, we conclude that, for (X, Y, Z) G M e , (VF(I, Y, Z), n) = 0. □ 

Let L°(X), for fc = 1,2,3, be the linear functions satisfying L^(0) = ifc, ^(1) = fk for fc = 1, 2, 
L^O) = 4>q and £3(1) = 0. Lemma [3.11 guarantees the existence of functions L k (X, Y, Z, e) for 
k = 1,2,3 such that for each e > and Y 2 + Z 2 < g 2 {X,e), L k (X,0,0,e) = L° k , L k (0,Y,Z,e) = 
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L° k (0), L k (l,Y,Z,e) = L° k (l), and J^L k (X,Y, Z,e) = when (X,Y,Z) £ M e . For each e > and 
k = 1,2, 3, introduce the functions L e k in terms of variables x, y and z: 

L U x ^y^ z ) = L k (X,Y,Z,e) = Lk(x,g(x,€)y,g(x,e)z,e). (3.2) 

Set 

u(x,y,z) =L\(x,y,z) - c x {x,y,z), 
v(x, y, z) =L\{x, y, z) - c 2 {x, y, z), 
4>{x, y, z) =L\{x, y, z) - $(x, y, z). 

Then, problem (|2.4p - (|2,5p is transformed into 

• {g 2 JJ T V(0 - L|)) = -A(ai(« - L\) - a 2 (v - L\)), 
^ = • {9 2 JJ T V(u - L\) - ai(u - Ll(x)) 5 2 JJW(0 - L§)) , (3.3) 

Yt =I gT V - {9 2 JJ T V{v - LI) +a 2 (v- Ll{x))g 2 jrV{cl> - L\)) , 
with the homogeneous boundary conditions: 

4>\lur =u\lur = v\lur = 0, 

3.4 

JJV) | M =(V«, JJV) | M = (Vv, JJV) | M = 0. 
System (|3.3p is supplemented with the initial conditions: 

u(0) = n , u(0) = v . (3.5) 

Introduce the subspace Hjj(Q) of H 1 ^): 

fll,(n) = {« € H\n) : U | LUR = 0}. 

Let M be the constant in Proposition 12.21 and let E e be the subset of Hj-,^) x Hp(Q) given by 

S e = {(«, u) e #|>(0) x fl^(n) : ai-Lf - M < a\U < a x L\, a 2 L% - M < a 2 v < a 2 L e 2 }. (3.6) 

It follows from Proposition 12.21 that, if (uq,vo) £ S e , then (u(t), v(t)) £ S e for every t > 0. 
Throughout this paper, for every e > 0, we denote by S e (t)t>o the solution operator associated 
with problem (|3.3|) - (|3.5p . We will use the same symbol A e to denote the global attractors of 
S e (t)t>o and that of problem (|2.4p - (|2.5p when no confusion arises. 
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The corresponding one-dimensional limiting system (|2.6p is transformed into 




(3.7) 



with the homogeneous Dirichlet boundary conditions 



<£, = u = v = 0, x = 0, 1 



(3.8) 



and the initial conditions 



u(0) = «o, and v(0) = vq. 



(3.9) 



Since So is an invariant region for problem (|2.6p - (|2.7p . we find that the one-dimensional problem 
(|3.7p - (|3.9p also has a positively invariant region which is given by 

E = {(u,v) G#o(0,l) xH£(0,1) : a x L\ - M < am < aiL?, a 2 L Q 2 - M < a 2 v < a 2 L° 2 } . (3.10) 

Similar to system f)2.6|) - (|2.7j) . problem ()3.T[) - (|3.9j) is well-posed in So, that is, for each (uo,t>o) £ So, 
there exists a unique solution (u, v) for problem (|3.7p - (|3.9p which is defined for all t > and (u, v) G 
C([0, oo), So). Further, the solutions are continuous in initial data with respect to the topology of 
-£fo(0, 1) x Hq(0, 1). Therefore, there is a continuous dynamical system S°(t)t>o associated with 
problem (|3.7p - (|3.9p such that for each t > and (uq,vq) G So, S°(t)(uo, vq) = (u(t),v(t)), the 
solution of problem (|3.7p - (|3.9p . When no confusion arises, we use the same symbol Aq to denote 
the global attractors of S°(t) t > and problem ff2T6|) - . 

3.2 Uniform estimates of global attractors 

In this section, we derive uniform estimates of the global attractors A t in e which are necessary 
for establishing the upper semi-continuity of A t at e = 0. In what follows, we reformulate problem 
(|3.3|) - (|3.5|) as an abstract differential equation in {{^(Q) x £[^(£1). 

Given e > 0, define an inner product (•, -)# e on L 2 (Q) by 




—^vw dx dy dz, 
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In the sequel, we denote \\w\\ p the standard norm of w for w G L p (0) orraG L p ([0, 1]), the 
standard norm of w for w G H S (Q) oimj£ -£P([0, 1]). Also, denote if e the space L 2 (J7) with the 
inner product (•, -)jj e , and X e the space Hj-)(Q) with the norm 



1 1 ^ 1/2 

= ( ||Vw|| 2 + ^ Ikylli + ^ll^lll 



Since Poincare inequality holds in Hj-,(Q), the norm ||u>||x E for w G Hjj(Q) is equivalent to the 
norm given by 

, 1.2 1„ „2 1„ l|2 XV2 

Miff* +^IKIl2 + ^IKIl2 

Due to assumption (|2.ip . there exist positive constants C±,C2, C3 (independent of e) and ei such 
that, for all < e < ei and a; G (0, 1), 

— < Ci, C 2 < ^ < C 3 . (3.11) 
5 e 



Consequently, y a e (u;, is equivalent to the norm ||u;||x £ , that is, 

C4wf Xe < a £ (w,w) < C 5 \\w\\ 2 Xe (3.12) 

for some constants C4 and C5 (independent of e). It follows from (|3.12p that for each e > 0, the 
triple {Hpifl), H e , o e (-, •)} defines a unique unbounded operator C e on Hj-,(Q) with domain D(C t ) 
in the following way: an element u G {{^(Q) belongs to D(C t ) if a e (w, u) is continuous in v G Hp(Q) 
for the topology induced from # e and (C € u,v)h £ = a e (u,v) for (u, v) G D(C € ) x Hj^^Q). In fact, 

L>(£ e ) = {u G #1,(0) : £ £ m G H e }, 

and for every u G D(C e ), 

C e u = -^rV • (g 2 JJ T Vu) . 

1/2 

Since the operator C e is self-adjoint on H e and positive, the fractional power C e is well-defined 
with domain D{C\ /2 ) = H^(Q), and for u G #2,(0), 

1 

H^I^Hik = a t (u,u). 
In view of f)3. 12|) there exist Cq and C7 such that 

OslMk < llAlk < C 7 ||n|U £ - (3.13) 
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With the above notations, system (|3.3p can be rewritten as 

C e <j> = XaAu - L\) - Xa 2 (v - L%) - 4rV • (g 2 JJ T VL\) , 

9 

+ DxCeU = -^-V • (g 2 JJ T VL\ + a x {u- L\)g 2 JJ T V(<j) - L%)) , (3.14) 
^ + D 2 C e v = -^V • {g 2 JJ T VL\ - a 2 (v - L\)g 2 jrv(<\> - L§)) . 

By the construction of functions L e k (k = 1, 2, 3), there exists ei > such that for any < e < e\, 
the following uniform bounds in e hold: 

Halloo + ||i|lk + \\J T VLl\\ Ht + || JJ T VL%\\ He + ||lv ■ {g 2 JJ T VL%) \\ He < C, (3.15) 

9 

where C is independent of e. Then it follows from the positive invariance of E e that there exists a 
constant C (independent of e) such that for any initial datum (uq,vq) € S e , the solution (u, v) of 
problem (|3.3p - (|3.5p satisfies, for all t > 0: 

IKt)l|oo + IKt)||oo<C and \\u(t)\\ He + \\v(t)\\ He <C. (3.16) 

Next, we start to derive uniform estimates of solutions in e in the space Hjj(Q) x Hjj(Q). 

Lemma 3.2. There exist a constant C (independent of e) and e\ > such that for any < e < ei 
and (uq,vo) E T, t , the solution (u,v) of problem 113. ,3\) - [3~5\) satisfies, for all t > 0: 

t+i 

(\Ht)\\ Xt + \\v(t)\\x.)dt<C. 
Proof. Taking the inner product of the first equation in (|3.14|) with <j> in H e , we find that 

\\CUf He = X ai {u-Ll4>) Hc -\a 2 {v-L\,<t>) Ht + (J T VLl,J T X7<P) He . 
By (fBTToT) and (I3T6D we have 

\\Cn\\ 2 He <Aai(||u|| He + Wi\\hM4>\\h 6 + Xa 2 (\\v\\ He + \\L\\\ He )\\4>\\H t + || I%\\ Ht \\cU\\H t 
<C\\cU\\H t <\\\ch\\ 2 Ht + \c 2 , 

which implies that 

\\CU\\h s < C. (3.17) 
Now, taking the inner product of the second equation in (|3.14p with u in H e , we get 

l^lHk + Di\\clu\\ 2 Ha = D x (J T VL\, J T Vu) H( + D iai ((« - L\)J T V{<i) - L|), J T Vu) He . 
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It follows from (|3.15p - (|3.17p that the right-hand side of the above is bounded by 

Ci||J T VLf|| ff J| J C|«|| ff ,+Diai(||«|| 00 + ||L!IU)(||^?^lk + II^VLIIkJH^IullH, 

<c\\clu\\ He < ^wchw^+d. 



Therefore, 



Similarly, 



±\\ u f Ht+ D 1 \\clu\\ 2 He <C 2 . (3.18) 
j t \\v\\ 2 He +D 2 \\clv\\ 2 He <C 3 . (3.19) 



Hence, for all f. > 0: 



| {hf He + IMiy + c * (\\clu\\ 2 Ht + \\jolvf H ^j <c 2 + c 3 , 

which, along (|3.13p and (|3.16p . implies Lemma 13.21 □ 

Lemma 3.3. There exist positive constants e\ and C such that for any < e < e\ and (uq, vq) 6 £ e , 
the solution (u,v) of problem \3.3\) - (375\) satisfies, for all t > 1: 

\\c e m\x. + \Ht)\\ Xt + Ht)\\x.<c. 

Proof. By (|3.15p . (|3.16p and the first equation in (|3.14p we get 

\\£e<f>\\H. < C (\\u\\ He + \\v\\ He + \\Ll\\ He + \\L 2 \\ H( + ||^V • ( ff 2 JJ T VL|) \\ He ) < C. (3.20) 
Taking the inner product of the second equation in (|3.14p with C e u in H e , we find 



(E±?^ r< ,..2 , rr 

V 9 2 



V-((u-L{)g 2 JJ T V(c/>-Ll)),C e u) . (3.21) 



By (]3.15p . the first term on the right-hand side of (|3.2ip is bounded by 

D l«-y , .2 t TTV7re\ r „, 1 I ^ n. II 1 V7 t„1 T TTX7Te\\\ . .11 r „.ll„ ^ 1 n. II r „.l|2 



^-V-(^JJ r VL|),£ eU I <Z) 1 [|- 5 V-(^JJ T VLf)|k.||A«||H e < -DiH^II^+C. (3.22) 



For the second term on the right-hand side of (|3.2ip . we have 

= -D iai (V(ti - • JJ T V(0 - Ll),Ceu) He 



D iai ( (u-L^V-( 5 2 JJ T V(0-^)),£ e u) . (3.23) 
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Using (|3.15p and (|3.20p . the first term on the right-hand side of (|3.23[) is bounded by 
D iai \ (V(u - L\) ■ JJ T X7(<P - L e 3 ),C e u) He \ 

< Diqi||V(u - L{) || 3 1| |r JJ r V(</> - L|)|| 6 ||£ e u|| 2 

< C\\V(u - L\)\\\ \\V{u - Ll)f Hl \\^JJ T V(0 - L\)\\ m \\C t u\\2 

< (\\Clu\\ Ht + \\J T VL\\\ H )j 2 [\\C e u\\ He + WFVLIWh. + II^V • (g 2 JJ T VL\)\\ H ) * 

x (\\CMH t + \\jrvLi\\ Ht + ll^v • (g 2 jrvLi)\\ H \ \\C e u\\ He 

1 
_> 



< C [\\C!u\\ Ht + CJ (\\C € u\\ He + C) 5 \\CM\h. 
^D^uf^ + CWcUw^ + C. (3.24) 
The second term on the right-hand side of (|3.23p can be estimated as 

1 



(u - L\)-V ■ (g 2 JJ T V(<P - LI)), C e u 

9 J H f _ 



< D iai (IMloo + IlLflloo) II^V • {g 2 JJ T V{4> - L%))\\ Ht \\CM\H t 



< D x a x (||«||oo + Halloo) (\\CM Hl + II^V • {g 2 JJ T VLl)\\ H }j \\C e u\\H e 

< C\\C e u\\ Ht < ~I>i||£ e u||ff« + C. (3.25) 

o 

Combining the estimates (|3.23p - (|3.25p . we obtain 

'D 1 Otl„ „ Tf , 2 t TW , , r *ss n \ , „ 1 

H, 



V ■{{u-Ll)g 2 JJ T V{<t>-Ll)),C e u) \ < -D 1 \\C e u\\ 2 Hi + C\\C}u\\ 2 H ^ + C. (3.26) 



It follows from (I3T2TT) . (I3T221) and (f3T26]) that, for all i > 0, 



Similarly, for all t > 0, 



|||>C!«||L +£> 1 ||£^||l f6 < CillAll^ +C 2 . (3.27) 



^||A|& e +£ 2 ||A«|&« <Ci||£?«|& e + C 2 . (3-28) 



Hence, we have, for all t > 0, 

||£ju|& e + ||£f Hlf^ + C 3 {\\C e uf He + \\C e vf He ) < Ci ( IIA^H^ + ||£!HI^ ) + C 2 , (3.29) 
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which, along with Lemma 13.21 and the uniform Gronwall's lemma, implies that, for all t >1, 

\\/:lu(t)f He + \\/:lv(t)\\ 2 He <c. 

The above estimate and the first equation in (|3.14p conclude the proof. □ 

Applying Gronwall's lemma to (|3.29[) for t E (0, 1), then by Lemma 13.31 and the first equation in 
(|3.14p we find that there exists e\ > such that, for any R > 0, there exists K depending on R 
such that for any < e < t\ and (uo^o) 6 S e with ||(wO)^o)||x e xX e < R, the following holds: 

\\C e (j){t)\\x t + IK*) Ik + IK*)lk < K, for t > 0. (3.30) 

An immediate consequence of Lemma [3. 3 1 also shows that all the global attractors A e are uniformly 
bounded in e in the space Hj-)(Q) x Hj^^), that is, the following statement is true. 

Proposition 3.4. There exist positive constants e± and C such that for all < e < e\ and (u, v) E 
A e , the following holds: 

\\(u,v)\\ XeXXc <C. 

The following is an analogue of Lemma 13.31 for the limiting system (|3.7p - (j3.9p . 

Lemma 3.5. There exists C > such that for any (uq,vq) £ So, the solution (u,v) of problem 
[iny-EBD satisfies, for all t>\: 



\\u(t)\\ m + \\v{t)\\ m <C 
In addition, there exists K depending on R when \\(uq,vq)\\ui x hi < R such that for all t > 0: 

\\u{t)\\ m + \\v{t)\\ m < K. 

Next, we establish estimates on time derivatives of solutions for both the three-dimensional 
system and the one-dimensional limiting system. 

Lemma 3.6. There exists ei > such that for any R > 0, there exists K depending only on R 
such that for any < e < ei and {uq,vq) G S e with ||(itO)^o)||x E xX E < R, the solution (u,v) of 
problem 113. 3\) - l[375\) satisfies 

pfur^u* o-ii^H 2 x H^H 2 ^ -i- 2 Alalia , n dv \\2 \ . < K Kt , >n 
Proof. Denote by 

~ d<p du dv 
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Differentiating (|3.14p with respect to t, we get 
C e (j) = Xa\u — Aa 2 £>, 



— + £>i£ e u = --^V • ^iVJJ r V((^ - L|) + ai(n - L\)g z JTV^ , 
+ £> 2 ^ = ~|V • (a 2 % 2 JJ r V(0 - Lf) + a 2 (v - L\)g 2 JJ T Vpj . 

From the above system, one derives 

C e (t4>) = \a\tu — \a2tv, 
— (in) + DxC e (tu) = u- -^V • (ai(tn) 5 2 JJ r V(0 - L|) + ai(u - L\)g 2 JJ T V(t$) 



j^(iS) + L> 2 £ e (i5) = v + • (a 2 (tv)g 2 JJ r V(0 - L|) + a 2 (« - L 2 ) 5 2 JJ T V(t0)) . 
The first equation in (|3.3ip gives 

\\C £ m\ He <C(\\tu\\ He + \\tv\\ He ). (3.32) 
Taking the inner product of the second equation in (|3.31j) with tu in H e , we have 

\f t \M% € +D l \\cl(tu)\\% e =D 1 ai J ^turv(<p-Li)-rv(tu) 

+ Dxai J ^{u - Ll)J T V{t4>) ■ J T V(tu) + t\\u\\ 2 He . 
By (|3.30p . the first term on the right-hand side of (|3.33[) is bounded by 
C|H| 3 || J T V(0 - L|)|| 6 || J T V(tu)\\ 2 < C\\tu\\l \\tuWl, ||J^V(0 - L§)[| Hl ||.TV(tu)|| 2 

< C\\tu^\\Cl{tu)\\l e {\\CMh. + II^IIIh.) (3-34) 

< c\\tu\^\\cl{tu)\\l e < \D x \\a{tu)f Ht + c\\mf HK . 



By ()3.32p . the second term on the right-hand side of (|3.33p is less than 

c\\u - L\\u\rv(t4>)\\ 2 \\rv(tu)\\ 2 < ^wchmk + c\\cl(t$)f Ht 

<^||4(«)[||r. + C'(||«[||r. + INlSr.)- 



(3.33) 



(3.35) 



Multiplying the second equation in (|3.14p by tu, after simple computations, we find that the last 
term on the right-hand side of (]3.33p satisfies 



t\\u\\ 2 Hf <C\\Cl{tu)\\ Hc { \\Clu\\ He + \\ClL\\\ He + ||£|<fc + \\CIL\\\ H€ 

(3.36) 

<l Dl \\cl(tu)\\ 2 Ht +C. 
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Combining the estimates in (|3.33p - (|3.36p . we get 

j t \\t^+D l \\Cl{tu)\\ 2 He < C(\\tuf He + \\tvf He )+C. (3.37) 

Similarly, 

j t \\tvf He +D 1 \\cl(tv)f He <C[\\tu\\ 2 He + \\tvf He )+C. (3.38) 
Finally, from ([3737D - (1338]) . we have 

j t {\\tu\\l e + \\tvf He ) + Ci (\\cktu)\\l, + ll^(^)ll^) < C (\\tu\\ 2 He + \\tv\\ 2 Hc ) + C, 
which, along with Gronwall's lemma, concludes the proof. □ 

We now describe the analogue of Lemma 13.61 for the one-dimensional limiting system (|3.7p - (|3.9p . 

Lemma 3.7. Given R > 0, there exists K depending only on R such that for any {uq,vq) 6 So 
with || («o, vq)\\hi x iji < R, the solution (u,v) of problem \3. 7| )-( TO|) satisfies 

MII w IIh 2 + II^II 2 + II^II 2 ) + / Mll^ll^ + ll^ll^)^<^ , t>o. 





Proof. The proof is similar to that of Lemma 13.61 but simpler, and therefore omitted here. □ 
3.3 Upper Semicontinuity 

In this section, we establish the upper semicontinuity of global attractors A e at e = 0. We first 
compare the solutions of the three-dimensional problem f|3.3[) - ()3. 5[) and the one-dimensional limiting 
problem (|3.7p - (|3.9p . and then establish the relationships between the global attractors of the two 
dynamical systems. 

In what follows, we reformulate limiting system (|3.7p as an operator equation. Let Hq be the 
L 2 (0, 1) space with the inner product (•, -)# given by 

{u,v)h = / gluvdx, 
Jo 

and let ao(-, •) be the bilinear form on (i^o(0, l)) 2 : 

. . / dw\ du)2\ f 1 o dwi dw2 , 

a Q {wi,w 2 ) = I -j— ,-r- =/ -7—dx. 

\ dx ax J Ho Jq dx dx 

For / G L 2 (Q), let M(f) G L 2 (0, 1) be the function: 



(M(f))(x) = - I f(x,y,z)dydz. 



7T 
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Lemma 3.8. Suppose f E H 1 ^). Then we have 

\\f-M(f)\\ He <Ce\\f\\x £ . 



Proof. Notice that 



l2 < x r , i 



f(x,y,z) / f(x,u,v)dudv 

7T 



2 



Using the identity 

/■</> ^ 

f(x,rcos9,rsin9)=f(x,pcos(j),psm(j)) — / —f(x,pcost,psmt)dt 

Je ut 

— / — — fix, t cos 0, t sin 9)dT, 

Jr OT 



dydzdx. (3.39) 



one can write 

-1 /-27T 



~ I I f (x,r cos 9, r sin 9)rdrd9 = f(x, p cos (f>, p sine/)) 
vr J r =o Je=o 

1 Z" 1 /" 27r / [4 ( df df \ \ 

/ / / — p sin t—(x,p cost, p sin t) + p cost— (x,p cost, p sin t) ) dt ) rdrd9 

7T J r=0 J e=0 \J t =e \ oy dz J J 



cos9—(x,Tcos9,Tsm9) + sin#— (x, r cos 9, r sin#) 1 dr ) rdrd9 



i r r w / 

7T 7 r=0 7^=0 \/r= 



<9y cfe 
Then, after simple computations, Lemma 13.81 follows from (|3.39p and the above. □ 

Let (?/>, P, Q) E (Hjj(Q)) 3 and let ((j) e ,u e ,v e ) be a solution of system ( |3.3I) . In view of the boundary 
condition (|3.4ft and the choices of L e k for k = 1,2, 3, we have 

-a e (0 e - = - \ai(u t - L\,ip) Ht + \a 2 {v e - Ll,ip) Hc , 

157 (l^' P ) H . = " °« ( "« - L! ' P) + (<* - LD£ ' /2( * - i5) ' £ ' /2P ) tt • (3.40) 

~k & ,Q ) H = _ <1,<t '' _ L '' 0) _ 0,2 ( ( "' _ i5)£ ' /2(A _ L ' ) " c ' 1/2<3 )„, ■ 

Let {4>,u,v) be the solution of the limiting system (|3.7p . View (<j),u,v) as an element in (-f/^fi)) 3 . 
Then a direct computation yields that, for (ip,P,Q) E (-ff^(r2)) 3 , 

-o 6 (^-L§,V)=-Aai(u-L? J V)ff. + Aa 2 («-ii§ J ^K+F(0-L§,V), 

^ p) = - a e (u - Ll P) + aa ((« - Li)Cy 2 ((f) - L^C^P)^ 

+ G 1 (u-Ll,<f>-LlP), (3.41) 

+ G 2 (v-L° 2 ,4>-L° 3 ,Q), 
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where, for appropriate functions p, q and r, and for i = 1, 2, 
F iP, Q) = [[ ~2 — )Px,q) + —Px,yq y + zq z 



Gi(p,q,r) = - ( ( - Px ,r) - ( 9 ^p x ,yr y + zr z ) (3.42) 

V V 9 % J J h € V 9 J H t 

+ { - l)l+lai {{^-^gf) Pq ^ r ) H +(-!) m ^ (jP^yry + zr^j H . 

Let 

= ^ - L\ - (0 - L°), P £ = Ue - L e x - (u - L°), Q e = v e - L\ - (v - L° 2 ). (3.43) 

Upon subtracting (13.41R from (|3.40p . we obtain that for any (■(/>, P,Q) G iP^,(f2) 3 , 

a e (V e ,V) =Aa 1 (P e ,^k e - Aa 2 (Q e ,^k e +F(0-L°,^), (3.44) 
i- (9 t P e , P) Hc = - a £ (P e , P) + ai ((u e - L\)cjip e , £f p) 

+ ai (V/;J(0-L°),£fp) -Gi(n-L?,0-L°,P), (3.45) 

^ (d t Q e ,Q) Ht = - a £ (Q £ ,Q) - a 2 ((u e - L 2 )£ e V,£|Q^ 

-a 2 (Q e Cl{4>-Lt),ClQ^ -G 2 {v-Ll,4>-Ll,Q). (3.46) 

For the above system, we have the following estimates. 

Lemma 3.9. There exists e\ > such that, for any R > 0, there exists a constant K depending on 
R such that, for any < e < ei and (uo,t>o) £ wi/i ||(uo, vo)||x e xX e < R, the following holds: 

\\P e (t)\\ 2 He + \\Q e (t)f He + W(t)\\h + f (\\P e (s)\\ 2 Xc + \\Q'(s)\\ 2 Xt ) ds < eKe K \ t > 0, 

Jo 

where (cj) e ,u t ,v e ) is the solution of problem ll3.3\) - [3~5\) with the initial condition (uo,vo), (4>,u,v) 
is the solution of problem \3. 7| )-( TOj) with the initial condition (M(uq), M(vq)), and (ip e ,P e ,Q € ) is 
given by j3.43\ ). 



21 



Proof. It follows from that 

a^,^) =Aai(P e ,^)// E -Aa 2 (Q e ,#)// £ + F(<f> - L° 3 ,if> e ), 



1 



(d t P*,P*) Hc =-a e (P*,P*)+ ai {u e -L\)£i^,CiP 



+ a l [P*a(4-L»),£lP e ) - Gi(u — L\, <j> - L 3 , P e ), 



(3.47) 



(3.48) 



D 2 



{d t Q e ,Q e ) Ht = - a e (Q e ,Q e ) - a 2 (« e - L e 2 )Ciijj e ,£iQ 



H, 



a 2 ( Q e £f (</> - L»),CIQ* ) - G 2 ( V - L°, - L°, Q e ). 



(3.49) 



Next, we estimate each term on the right-hand sides of (j3.47|) - f)3.49|) . The first two terms on the 
right-hand side of (|3.47|) are bounded by: 



\ ai |(P e , V £ kJ + Aa 2 |(Q e , #) He | <C(||P1h £ + IIQIk) ll^lk 



(3.50) 



By (|3.1ip we find that g satisfies 



d x g 2 d x gl 



„/,el 



< Ce. 



Then, by Lemma 13.51 the first term in F((f) — L®,tp € ) on the right-hand side of (|3.47p is less than 



d x g 2 d x gl 



9l 



He 



<Ce\m\^<Ce' + -a t {r^ e )- 



(3.51) 



It follows from (|3,30p and Lemma [3 . 5 1 that the second term in F((p — L 3 , tp e ) on the right-hand side 
of (ET4T)) is bounded by, for t > 0, 



^{cp - L° 3 ) x , yd y (ct> € - L|) + zd z (<f> e - LI) 
9 / H e 

< C (\\d y <f> e \\ He + R&lk + \\dyL e 3 \\ He + < Ce. 

By (jlLlZD and (|330>(|332l . we obtain, for all t > 0, 

\mt)\\h < C«e(#,#) < C (||P £ |& e + + Ce. 



(3.52) 



(3.53) 



We now deal with the right-hand side of (|3.48p . By (|3.53p . the second term on the right-hand side 
of (ET38]) is less than 



{u e -Lx)LW,£,lP l 



< C(W\\ 2 x, + W^Wk) < C(\\P%ic t + WQTx.) + Ce. (3.54) 



3>e||2 



ei|2 



ie||2 
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Since the functions <fi and L3 depend on x G (0, 1) only, we have 

cu = «/ T v</> = (<9^,o,o) T , CHI = J t ^lI = (a,L°,o,o) r , 

which, along with Lemma 13.51 and the first equation of (|3.7p . implies that, for all t > 0, 

Halloo = Halloo < c\\d x <f>\\ H1 < c\\<p\\ H 2 < c. 

By (|3.55p . the third term on the right-hand side of (|3.48[) is bounded by 



(3.55) 



< ai \\C 



\P 6 



< c\\n\^ + -a e (p^n. 



(3.56) 



Note that the term G\ on the right-hand side of (|3.48[) can be estimated in a similar manner as 
(|330 >(|332|) . Therefore, it follows from (f3^5D and ^M)-^M that, for t > 0, 



J t W p£ Wk + W pe Wx e < C(\\P e f He + \\Q% C ) + Ce. 



Similarly, Q e satisfies, for t > 0, 



d 



f t \\Q e \\k + \\Q e \\x. < c(\\n\ 2 He + iiQiy + ce. 



Then, it follows from ([3371) - (13381) that, for t > 0, 



fit 



P e \\k + Uplift) + H^llxe + \\Qlx c < C(\\P*f He + HQlU + Ce. 



By Gronwall's lemma, we get 

ll^Wlllr. + \\Q e (t)\\k < z ct (ll^(o)llk + \\Q € m 2 Hc ) +^ ct 

< Ce ct {\\u - M(u )f He + \\L\ - + \\v - M(v )\\ 2 Hc + - L°|&J + ee 

By U32J we see that L U L 2 G PF 1>0 °(J2 e ), and hence, for k = 1,2, 



ire r0||2 



yg^r(x, sgy, sgz) + z 9^-{x, sgy, sgz) ) ds||^ e < Ce 2 



From (I3.60l) - (l3.6ip and Lemma f3T8l we find that 

\\P e {t)f He + W{t)f He <e{C + l)e ct 
Integrating (|3.59p between and t, by (|3.62p we conclude Lemma [3? 



(3.57) 



(3.58) 



(3.59) 



(3.60) 

(3.61) 

(3.62) 
□ 
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Next, we improve the uniform estimates in e given in Lemma 13.91 

Lemma 3.10. There exists e\ > such that, for any R > 0, there exists a constant K depending 
on R such that, for any < e < ei and (wo,t>o) £ E e with ||(uo, Vo)||x e xX e < R, the following holds: 

,dP € n2 u dQ E 



dt 



\k + H^IIh.) + * (W pe Wx. + IIQ e ll*.) < v^**, t > 0, 



where (ip e ,P e ,Q e ) is given by [3.4-3fy , (<f> e ,u e ,v e ) is the solution of problem \3. 3\) - {375\) with the 
initial condition (uq,v ), and (cp,u,v) is the solution of problem \S. 7\ )- [3\V\) with the initial condition 
(M(u ),M(v )). 

Proof. Denote by 

P*-^- fr-^ fr-W- (3 63) 

F ~ dt ' Q ~ dt ' ^ ~ df {6 - b6) 

Differentiating systems (j3.44f l (|3.46p with respect to t, multiplying the resulting systems by t, 
replacing ip, P and Q by tip, tP and tQ, respectively, we obtain 

a £ {t*p e ,t4> e ) =\ ai (tP e ,t4> e ) Ht - Xa 2 {tQ e M e ) Ht +tF(4> t ,t-ip e ), (3.64) 



1 d 
2LhJt 



-^\\tP e f He +a e (tP e ,tP e ) =a x t (d t u e Ct^,c\tP e ^j +a x t (^(u e - L\)CJi) e , £f tP^j 



i \ /i 



+ ait ( P e Cl (0 - L° 3 ),C!tP e + ait P e £f^,£| tP e 

J H e V / H e 

ml I (^f-^f) (u-LD^tP^ (3.65) 



9 2 9o J J H e 



- ait [—(u - L»)(j) tx ,tyd y P £ + tzd z P e 

- tG x {ut, 4> - L°, tP") + -L(P%tP e ) He , 

^j t \\tQ e \\% e + a e (tQ\tQ e ) = -a 2 t (d t v e /:h* , £? tfrj - a 2 t ^(v t - U 2 )Cl^,CjtQ^ 

- a 2 t (Q e Cj{4> ~ L° 3 ),dtQ e ) - a 2 t (q* C} C}tQ e 

+ a 2 t ( (^4- ~ ^f) (« " L° 2 )<j> tx ,tQ<) (3.66) 
\\ 9 % J J He 



+ a 2 t — (v - L 2 ) )4> tx ,tyd y Q t + tzd z Q € 

V 9 J H ri 

- tG 2 (v t , - L° 3 , t&) + -L(Q e , t&) He . 

P>2 
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We now estimate every term involved in the above system. Note that (|3.64p implies that 

\\cltr\\ 2 He < C (\\tP e \\ 2 He + \\m\ 2 Hc ) + Ce 2 t 2 . (3.67) 
By (I3.30P and Lemma 13.91 we see that the first term on the right-hand side of (|3.65|) is bounded by 

\ ai t (dtuxIrxitP^ l < ct\\d t u e \\ 6 \\^h £ h\\£ltp% 

<C't||9 t « e ||^i||£!#|||||£!^||| 1 ||£!tP e || 2 

< ±\\cjtP% c + ct 2 \\cld t u e \\ 2 He \\cl^\\ He \\£^ e \\ He 

< ^\\dtP% e + V~eCe ct \\td d t u e \\ 2 He . (3.68) 
By (j3.67p . the second term on the right-hand side of (|3.65p is less than 

a x t {{u e - L\)cir,£}tP^j < ±\\dtP% ( + C (\\tP e \\ 2 H( + \\tQ £ f He ) + e 2 Ct 2 . (3.69) 
By Lemma l3.7( the fourth term on the right-hand side of ()3.65|) is bounded by 
C\\tCU t \U\P e \\ 2 \\cltP% < L\\£,ltP e \\ 2 Ht + C\\t<f> t \\ 2 H 4P e \\H e < ^\\chn\ 2 Ht +eCe ct . (3.70) 

Other terms on the right-hand side of (|3.65p can be estimated in a similar way as the proof of 
Lemma [3791 Therefore, by (|3.65p . (|3.68p - (|3.70p and the estimates for other terms, we have 

a^H^Hir. + \\\^ tpe \\H, < C (\\tP e \\h + WtQlk) + zCe ct + e 2 C {\\tu t \\ 2 m + \\tv t f m ) 

+ V~eCe ct \\t£ld t u e \\ 2 He + -L(P*,tP e ) He . (3.71) 

Next, we deal with the last term on the right-hand side of the above inequality. Replacing P in 
([SagD by td t P e = tP € , we get 



l-(pt,tpt) He + l t j t \\CjP e \\% = ai ( (», - L) )£,V . C?IP' 



+ ai (^P e cl(<P-Ll),citP e j -G x {u-L\,(j)-Ll,tP e ). 
Using Lemma 13.91 and proceeding as before, we obtain from the above that 

l-^,tn^ + ^tj t \\cln\ 2 He < \\\chn\ 2 Hf +eCe ct + e 2 C. (3.72) 
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Then it follows from ([3T?ll - ([3T72l ) that 



1 d 1 1 1 d 1 

2L>idt" 2" 2 dt" UHe 

<c(\\tn\^ + \\tQ^)+eCe ct 

+ £ 2 C {\\tu t \\ 2 m + \\tv t f m ) + V~eCe ct \\tCldtu e \\ 2 He , 

which implies that 

\jt (^II^IIh. +*II^P e [|!r.) < \0P e f H , + c(\\tP*t He + +eCe™ 

+ e 2 C {\\tu t \\ 2 Hl + IKI^i) + V~eCe ct \\t£ld t u e \\ 2 He . (3.73) 
Similarly, by equation (|3.66p . we can show that 

\jt {j5$ tW »< +*H£?Q e |lk) ^ ^IQI^ + + ||tQl| e ) +eCe Ct 

+ e 2 C {\\tu t \\ 2 m + \\tv t \\ 2 H1 ) + ^Ce ct \\tdd t v e \\ 2 He . (3.74) 

By (I3T73D - (I3T71D we find that 

| ^\\tPi 2 He +t\\cin\ 2 Hc + ^pQ e nk +^iQ e n^ 

< c (±^tf*f Bt + t\\cin\ 2 Ht + ^- 2 \wf He + twdQlk) + \\£ip € \\k + W^Qlk 

+ eCe ct + e 2 C {\\tu t \\ 2 Hl + \\tv t \\ 2 Hl ) + JlCe ct (\\tcl d t u £ f He + \\tcl d t v e f H ) , 

which, along with GronwalFs lemma and Lemmas 13.61 13.71 and 13.91 implies Lemma 13.101 □ 

Let (c\ , c|, <J? e ) be the solutions of problem (|2.4p - (|2.5p with initial datum (ci,o, £2,0)1 and (ci, C2, $) 
be the solutions of problem (|2.6|) - (|2.7p with initial datum (M(ci i o), -^(c2,o))- Then as an immediate 
consequence of Lemma 13.101 we find the following estimates which are essential to prove the upper 
semi-continuity of the global attractors. 

Lemma 3.11. There exists e\ > such that, for any R > 0, there exists a constant K depending 
on R such that, for any < e < t\ and (01,0,62,0) S £ with ||(ci ( q, C2 ) o)||x e xX £ < R> the following 
holds: 

{Hi*) - Cl (t)\\ 2 Xe + \\c\{t) - c 2 (t)\\x e ) < VeKe Kt , t > 1. 
We are now in a position to prove the upper semi-continuity of global attractors. 
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Proof of Theorem 12.31 Let T e (t)t>o and T°(i)t>o be the solution operators of problem (|2.4p - 
(|2.5p and problem (12.6p - (12.7p . respectively. Then it follows from Proposition 13.41 that there is a 
constant R > (independent of e) such that 

||(ci,c 2 )||x £ xX £ < R, for all (c 1? c 2 ) G -A £ - 

For the given r\ > 0, since *4o is the global attractor of T°(t), there exists ro = tq(t],R) > 1 such 
that, for any t > tq, 

inf \\T (t)(Mz)-z \\x e xX e < 2 
for any z = (ci, 02) G ,4 e . On the other hand, by Lemma 13.111 we find that 

\\T*(t )z-T (t )(Mz)\\ xm < e*K{R)e K ^\ 

for some constant K{R). Therefore, we obtain that, for any z = (ci, C2) G ^4 e : 

inf ||T e (r )z - z ||x e xx £ < ? + ei^fl)^, 

which implies that, for e > small enough: 

dist XeXXe (T e (T )A e ,A ) < 77. 

The proof is completed since T^ro).^ = ^4 e . 

4 Steady-states for the one-dimensional limiting PNP system 

Many mathematical works have been done on the existence, uniqueness and qualitative properties of 
boundary value problems even for high dimensional systems and algorithms have been developed 
toward numerical approximations (see, e.g. |15^ \16[ [25l I17j). Under the assumption that fi <C 
1, the problem can be viewed as a singularly perturbed system. Typical solutions of singularly 
perturbed systems exhibit different time scales; for example, boundary and internal layers (inner 
solutions) evolve at fast pace and regular layers (outer solutions) vary slowly. For the boundary 
value problem (|2.8|) and (|2.9|) . there are two boundary layers one at each end. Physically, near 
boundaries x = and x = 1, the potential function <p{x) and the concentration functions c\{x) 
and C2(x) exhibit a large gradient or a sharp change. In [5], for a\ = a-i = 1, the boundary value 
problem for the direct (with h(x) = 1 in (|2,8p ) one-dimensional PNP system was studied using 
the method of matched asymptotic expansions as well as numerical simulations, which provide 
a good quantitative and qualitative understanding of the problem. In [23], geometric singular 
perturbation theory (see, e.g. [H [181 EH E2]) was applied to the study of this singular boundary 
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value problem. The treatment for the limiting one-dimensional PNP system (|2.8p carrying the 
geometric information of the three-dimensional channel follows that in [23] . 

It is convenient to study an equivalent connecting problem to the boundary value problem (|2,8p 
and (|2.9|) . Let Bl and Br be the subsets of R 7 defined, respectively, by 

Bl ={<P = <Po, v = -h(0)(aih - a 2 / 2 ), w = a\l\ + a\h, r = 0}, 

(4.1) 

Br ={</> = 0, v = —h(l)(ctiri - a 2 r 2 ), w = a\r\ + a 2 r 2 , r = 1}. 

The boundary value problem is then equivalent to the following connecting problem: finding a 
solution of fl2~T0D from B L to B R . 

For n > 0, let M£ be the union of all forward orbits of (12. 10h starting from Bl and let be the 
union of all backward orbits starting from Br. To obtain the existence and (local) uniqueness of a 
solution for the connecting problem, it thus suffices to show and intersect transversally. The 
intersection is exactly the orbit of a solution of the boundary value problem, and the transversality 
implies the local uniqueness. The strategy is to obtain a singular orbit and track the evolution of 
and Mr along the singular orbit. As discussed in the introduction, a singular orbit will be a 
union of orbits of subsystems of (|2.10p with different time scales. 

The boundary layers will be two orbits of f)2 . 14|) : one from Bl to Zq in forward time along the 
stable manifold of Zq and the other from Br to Zq in backward time along the unstable manifold 
of Zq. The two boundary layers will be connected by a regular layer on Zq, which is an orbit of a 
limiting system of (|2.10j) . The next two subsections are devoted to the study of boundary layers 
and regular layers. 

4.1 Fast dynamics and boundary layers 

We start with the study of boundary layers governed by system (|2.14p . This system has many 
invariant structures that are useful for characterizing the global dynamics. 

The slow manifold Zq = {u = v = 0} consisting of entirely equilibria of system (|2.14p is a 
5-dimensional manifold of the phase space M. 7 . For each equilibrium z = (4>, 0, 0, w, J\, J 2 , r) G Zq, 
the linearization of system (|2.14p has five zero eigenvalues corresponding to the dimension of Zq, 
and two eigenvalues in directions normal to Zq. The latter two eigenvalues and their associated 
eigenvectors are given by 

A± = ±v / ^ and n± = {{±^)~ l , 1, ±^fw, ±(a 2 - a^y/w, 0, 0, 0) r . (4.2) 

Thus, every equilibrium has a one-dimensional stable manifold and a one-dimensional unstable 
manifold. The global configurations of the stable and unstable manifolds will be needed for the 
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boundary layer behavior. For any constants J*, J 2 and r*, the set J\f = { J\ = J{, J 2 = J 2 , r = r*} 
is a 4-dimensional invariant subspace of the phase space M7 . 

Surprisingly, system (|2.14p possesses a complete set of integrals with which the dynamics can 
be fully analyzed; in particular, the stable and unstable manifolds can be characterized and the 
behavior of boundary layers can be described in detail. 



Proposition 4.1. (i) System {2.1$ has a complete set of six integrals given by 



«2 — ot\ a.ia.2 tt 1 m \oL\v/h(r) + w\ 
Hi = W — — v u , ±±2 — • 



h(r) 2/i 2 (r) ' z Y a 2 

In laov/hir) — w\ 

H 3 =<j) + 1 ' V ' ] -, H 4 = J h H 5 = J 2 and H 6 = r, 

ati 

where the argument of Hi's is ((f), u, v, w, Ji, J 2 , t). 

(ii) The stable and unstable manifolds W s (Zq) andW u (Zo) of Zq are characterized as follows : 

W S (Z ) = U{W s (z*) : z* G Z } and W U {Z ) = U{W u (z*) : z* G Z } 

and, for z* = 0, 0, w*, Jf, J|, r*) G Z , a point z = (<p,u,v,w,J 1 ,J 2 ,T) G W s {z*) U W u {z*) if 
and only if 

In w* In w* 

H 1 {z)=w*,H 2 {z) = <l>* ,H 3 (z) = cf>* + 



Q! 2 OL\ 
J\ =J*, J 2 = J 2 , T = T*. 

(Hi) The stable manifold W s {Zq) intersects Bl transversally at points with 



u = -[sgn (a 2 l2 - *,h)]S2h{J k + l 2 - ^ + ^(a.h)^ {a^)^ ^ 

V ai«2 

and arbitrary J\ and J 2 , where sgn is the sign function. The unstable manifold W u (Zo) intersects 
Br transversally at points with 



u = [sgn (a 2 ^2 - air 1 )\\/2h(l)\ n + r 2 (4.4) 

V aia 2 

and arbitrary J x and J 2 . Let N L = B L n W S (Z ) and N R = B R D W U (Z ). Then, 



to(N L 



{ ( 1 U\l\ Q 2 a l 

i </>o + 1 In — -,0,0, (ai + a 2 )(a 1 l 1 ) a i+ a 2 (a 2 Z 2 )<n+<*2 , h, ^2,0 

[ \ ct\ + a 2 a 2 / 2 



1 airi 



afJVfl) = <^ In ,0,0, (ai + a 2 )(airi) Q i+ Q 2 (a 2 r 2 ) a i+ a 2 , J 1; J 2 , 1 

a\ + a 2 a 2 r 2 



/or a// Ji and J 2 . 
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Proof. The statement (i) can be verified directly. The statement (ii) is a simple consequence of (i) 
together with the fact that 0(£) — ► <fi*, u>(£) — > w* , u(£) — > and v(£) — > as £ — > oo for the stable 
manifold and as £ — > — oo for the unstable manifold. 

For the statement (hi), we present only the proof regarding the intersection of W s {Zq) and Bl. 
Suppose 

z° = (<p°,u ,v ,w ,J?,J$,0) = (<Po,u ,h(0)(a2l2-a 1 l 1 ),alh+a 2 2 l2,J^Jl,0) 

is a point in Bl H PF s (-Zo). Then, using the integrals H\, H2 and H%, the solution z(£) = 
((/>(£), n(£), f(£), w(£), J°, J2, 0) of system f)2. 14|) with initial condition z(0) = z° satishes 

ffa(*(£)) =0(0 - — In |ai«(0//»(0) + u>(£)| = B, 

fl3(«(0) =<H0 + — In |a 2 «(e)//»(0) - = C 

ai 

for some constants A, 5 and C, and for all £. From the initial condition, we get 

ln(ai + a 2 ) ln(a 2 Z 2 ) , ~ , , ln(ai + a 2 ) , hi(aiZi) 
B = (po and G = (po H h 



a 2 a 2 ai ai 

Since u(£) — ► and u(£) — > as £ — > +00, we have that u>(+oo) = A and 

ai + a 2 ai + a 2 

C — B = in u;(+oo) = In A. 

«ia 2 aia 2 



Hence, 



Therefore, 



w(+oo) = A = (ai + a 2 )(ai/i) Q i+ a 2 (a 2 / 2 ) ai+a2 • 



r sgn (v°)]V2h(0)Jh +h~ — and 0(+oo) = O + In 



a\02 ai + a 2 «2^2 

The choice of the sign for u° comes from the consideration that the stable eigenvector n_ in (|4.2p has 
u and u components with opposite signs. Thus, Bl and W s (Zq) intersect at the points with u = u° 
given above, and all Ji and J 2 . IfiV L = B L DW S (Z Q ), then u(N L ) = {(<^(+oo), 0, 0, w(+oo), J x , J 2 , 0)}. 
The above formulas for 0(+oo) and u>(+oo) = ^4 gives the desired characterization of lo(Nl). Lastly, 
since the stable manifold is completely characterized, one can compute its tangent space at each 
intersection point to verify the transversality of the intersection. It is slightly complicated but 
straightforward. We will omit the detail here. □ 
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Part (iii) of this result implies that the boundary layer on the left end will be an orbit of (|2. 14[) 
from (4>o, ul, ct 2 l 2 — (Xih,a% l\ + ot\l 2 , Ji, 0) G Bl to the point 

(1 Ql/l "2 "1 \ 

00 + : In— -,0,0,(ai +a 2 )(ai/i) Q i+ Q 2(a2^) ai+Q2 ,«/i,«/2,0 £ Z , 
a\ + a 2 a 2 L 2 J 



where Ul is given by the display (14. 3D and I\ and I 2 are arbitrary at this moment; and that on the 
right end will be a backward orbit of (|2.14j) from the point (0, ur, a 2 r 2 — OL\r\, afri+a 2 r 2 , Ji, J 2 , 1) £ 
Br to the point 

(1 OL\T\ "2 "1 \ 
: In ,0,0,(0! + a 2 ) (am) a i+ a 2 (a 2 r 2 ) a ^ , Ji, J 2 ,l € Z , 
at + 02 a 2 r 2 J 

where ur is given by the display (|4.4p and Ji and J2 are arbitrary at this moment. It turns out 
that there is a unique pair of numbers J\ and J 2 so that the corresponding points zl and zr can 
be connected by a regular layer solution on Zq. The regular orbit together with the two boundary 
layer orbits provide the singular orbit. 

4.2 Slow dynamics and regular layers 

We now examine the slow flow in the vicinity of the slow manifold Zq = {u = v = 0} for regular 
layers. Note that system (|2.1'2|) resulting from (|2.1U|) by setting [x = reduces to u = v = and 

k = 0, J 2 = 0, t = 1. 

The information on cp and w is lost. This indicates that the slow flow in the vicinity of Zq is itself a 
singular perturbation problem. To see this, we zoom into an 0(^)-neighborhood of Zq by blowing 
up the u and v coordinates; that is, we make a scaling u = fj,p and v = fiq. System ()2.10j) becomes 

: 1 ■ h r( T ) f t 

<p =—■ -p, HP = Q, /J,q = pw + (J, — q + (a x Ji - a 2 J 2 ), 
h(r) h{T) 

a\OL 2 a 2 -a\ a\Ji + a 2 J 2 (4.5) 

ii =0, j 2 = 0, f = 1, 

which is indeed a singular perturbation problem. When ^ = 0, the system reduces to 

4>=-ri—rP, = q, = pw + (ai Ji - a 2 J2), 
/i(r) 

a 2 - ai a?Ji + ao<^2 (4.6) 

h(T) h(T) 

j\ =0, 3% = 0, f = 1. 
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Dynamics of <fi and w survives in this limiting process. For this system, the slow manifold is 



i oliJi — OL\J\ 

So = <i p = , q = 

uu 



(4- 



The corresponding fast system is 

$ = ^77~T> P =?> q' = pw + (aiJi - a 2 J2) + fJ--^~q, 
h{T) h{T) 

, o a i a 2 "2 - "i a?Ji + a?J 2 (4.7) 

M =M im PQ + ^^V PW " " Mr) ' 

J( =0, 4 = 0, r' = 0. 

The limiting system of (|4.7p when fi = is 

0' =0, p' = g, q' = pw + (ai Ji - 0^2), 
«/ =0, J[ = 0, J 2 = 0, r' = 0. 

The slow manifold So is the set of equilibria of (|4.8p . The eigenvalues normal to So are A±(p) = 
±^/w. In particular, the slow manifold So is normally hyperbolic, and hence, it persists for sys- 
tem (|4.7p for fj, > small (see [6]). 

The limiting slow dynamic on So is governed by system (14. 6h . which reads 

• a 2 J 2 -aiJi . «ia 2 (Ji + J 2 ) f n 1 

<P= rr~\ ' w = ZT~\ > Ji = °' T = l - 

n(T)w n{T) 

The general solution is characterized as: J\ and J 2 are arbitrary constants, and 

f x 1 

t(x) =to + x, w(x) =w - aia 2 (Ji + J 2 ) / — ds, 



h(r Q + s) 

X) =Vo + («2 J2 - OL\ Jl) / 7^ ; C?S (4.9) 



a 2 J 2 — ai Ji 

=zy ° / t I t \ In 

otiat2{Ji + J2) 



_ «i« 2 (Ji + J2) r 1 ds 

^0 Jo ^(jo + s) 



where to = r(0), (f)(0) = vo and w(0) = wq. Note that, if J\ + J 2 = 0, then = u>o and 

a 2 J 2 - aiJi P 1 
</>(£) = f o H / T7 : — ^ as- 



The latter is the limit of <j)(x) in (|4.9p as J\ + J 2 — > 0. We thus use the unified formula (14. 9j) even 
if Ji + J 2 = 0. 

To identify the slow portion of the singular orbit on So, we need to examine the cj-limit (resp. 
the a-limit) set of M£ n W S (S ) (resp. M% n W U («S )) as fi -» 0. To do this, we fix an 0(1)- 
neighborhood of Sq. In terms of U and V, this neighborhood is of order 0(/x). For /i > small, 
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the time taken in terms of £ for M^ and M^ to evolve to any 0(/i)-neighborhood of {u = v = 0} 
is of order 0(/j,\ ln//|). Thus, the A-Lemma ([5]) implies that M^ (resp. M^) is C 1 0(/x)-close to 
M£ (resp. M^) in any O (/^-neighborhood of {u = v = 0}. Therefore, in an 0(l)-neighborhood of 
So in terms of p and q, (resp. MjJ) intersects W s (5o) (resp. VF u (5o)) transversally. And, by 
abusing the notations, if N L = M£rW s (<S ) and N R = M^nW u (S ), then u(N L ) and a(iVfl) have 
the same descriptions as those in Proposition 14. II with u = v = replacing by p = {a.%Ji — a\J\)jw 
and q = 0. 

The slow orbit should be one given by (|4.9|) that connects uj(Nl) and a(Nji). Let M^ (rep. Mr) 
be the forward (resp. backward) image of cv(Nl) (resp. a(Nji)) under the slow flow (|4,6p . 



Proposition 4.2. M^ and Mr intersect transversally along the unique orbit given by ( f^.ff| ) /rom 
x = to x = 1 wdta 

r = 0, w = (qi + a 2 )(ai/i) a i +Q! 2 (a 2 ^) Ql+a2 , fo = 4>o H hi — -, 

and Ji and J 2 are as given in Theorem\2.4\ 



Proof. We show first that Ml and Mr intersect along the orbit with the above characterization. 
In view of (14. 9p and the descriptions for uj(Nl) and oi(Nr) in Proposition 14.11 the intersection is 
uniquely determined by 

a 2 "1 

t = 0, u>(0) = (ai + a 2 )(ai^i) a i +Q 2 (a 2 / 2 ) a i+«2 , 

"2 "1 

•u>(l) = (ai + a2)(a;iri) a! i+ a 2 (a2r 2 ) Q i+ a 2 f 
u = HO) = 00 + — ^- In 0(1) = — L- In 5™. 

Substituting into (|4.9p gives 



Oi\ -\- Oi2 ( a 2 a l "2 "1 \ 

di + d2 = 1 ( («i/i) a i+«2 (02/2) Q i +a 2 — (ain) a i+ a 2 (a^) a i+ a 2 J ; 

«i«2 Jo h~ 1 (x)dx \ ' 

( °>2 "1 "2 °1 \ 

(aih) a i +a 2 {oi2h) ai+c " 2 — (airi) Q i+ a 2 (a2r2) a i +a 2 J 

which in turn yields the expressions for J\ and d2- To see the transversality of the intersection, 
it suffices to show that u(Nl) • 1 (the image of uj(Nl) under the time one map of the flow of 
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system (|4.6j) ) is transversal to o(Nr) on So H {r = 1}. If we use (<j>,w, J%, J 2 ) as a coordinate 
system on So n {r = 1}, then the set uj(Nl) ■ 1 is given by {(0(Jx, J2), J2)) ^1, ^2)} with 

if-r j \ 1 1 i "1^1 a 2 J2-a\Ji , A Poaia 2 (Ji + J2) 

J2) =<f>o h ; — in — : rr~r~r\ m 1 



a\ + Ct 2 «2^2 ai«2(«/i + ^2) V 
w(Ji, J2) ={ol\ + a 2 )(ai/i) Q i+ Q 2 (a 2 l 2 ) a i+<*2 - Poa\a 2 (Ji + J 2 ), 

where po = Jq h~ 1 (x)dx. Thus, the tangent space to u)(Nl) ■ 1 restricted on So D {t = 1} is 
spanned by (0^,^,1,0) = (0j 13 -poa/3, 1, 0) and (0./ 2 , wj a , 0, 1) = (0j 2 , -poa(3, 0, 1). In view of 
the display in Proposition 14. H the tangent space to a(Nji) restricted on So n {r = 1} is spanned 
by (0,0, 1,0) and (0,0,0, 1). Note that So n {r = 1} is four dimensional. Thus, it suffices to show 
that the above four vectors are linearly independent, or equivalently, <pj 1 ^ (j)j 2 . The latter can be 
verified by a direct computation. Indeed, if J\ + J 2 / at the intersection points, then 

, , «i + «2 , ( P$oiia 2 {Ji + J 2 )\ 
<PJi ~ <Pj 2 = /j i j \ In 1 + 0; 



aia 2 {Ji + J 2 ) \ w 

if Ji + J 2 = at intersection points, then (j>(Ji, J 2 ) = vo+ po(& 2 J 2 — a\Ji)/wo and hence —cftj 2 = 
-po(ui + a 2 )/wo / 0. □ 



4.3 Proof of Theorem E31 

We provide a detailed version of Theorem 12.41 and its proof. 

Theorem 4.3. Assume that a\l\ ^ a 2 l 2 and ct\r\ ^ a 2 r 2 . For fj, > small, the connecting 
problem \2.10\) and \4- 1\ ) has a unique solution near a singular orbit. The singular orbit is the 
union of two fast orbits of system \2.1J$ and one slow orbit of system J^.6p ; more precisely, with 
both J\ and J 2 given in Theorem \2.4\ 

(i) the fast orbit representing the limiting boundary layer at x = lies on Bl n W s {Zo) from Bl 
to Lo{Ni) C Zq whose starting point has the u-component given by ^4-3^ in Propositions \4-l\ 

(ii) the fast orbit representing the limiting boundary layer at x = 1 lies on Br n W u (Zq) from 
Br to a(Nji) C Zq whose starting point has the u-component given by \4-4ty in Propositions \4-l\ 

(Hi) the slow orbit on So connecting the two boundary layers from x = to x = 1 is displayed 
in {4-ty together with the quantities in Proposition ^. 2\ 

Proof. The singular orbit has been studied in Sections 14.11 and 14.21 which is summarized in (i), (ii) 
and (iii) of this theorem. It remains to show the existence and uniqueness of a solution near the 
singular orbit for fi > 0. Recall that M£ (resp., M^) is the union of all forward (resp., backward) 
orbits starting from Bl (resp., Br). It suffices to show that, for \x > small, M£ and intersect 
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transversally with each other around the singular orbit. We note that the assumption a%h ^ (X2I2 
and a\r\ ^ «2 r 2 imply the vector field of (|2.10p is not tangent to and Br, and hence, and 
M R are smooth invariant manifolds. 

For (j, > small, the evolutions of M£ and M R from B^ and Br, respectively, to an /z- 
neighborhood of Zq along the two boundary layers are governed by system (|2,13p . Since, for 
system (pTilD . M° and M R intersect W S (Z ) and W U (Z ) transversally, we have that M£ and 
intersect W s (Zq) and W u (Zq) transversally. As discussed in Section B~2l in terms of the blow-up 
coordinates, and M R intersect W s (So) and W u (Sq) transversally for system (|4.7|) . And, if we 
denote N L = M^n W S (S ) and N R = M R n W U (S ), then the vector field on S is not tangent to 
lo(Nl) and o:(Nr). Furthermore, the traces Ml and M^j of uj(Nl) and o:(Nr) respectively under 
the slow flow on So intersect transversally. All conditions for the Exchange Lemma (see [30J and 
also [20l [18j CEH] ) are satisfied, and hence, and M R intersect transversally. The intersection has 
dimension 

dim M£ + dimAfg - 7 = 4 + 4- 7 = 1, 
which is the orbit of the unique solution for the connecting problem near the singular orbit. □ 

Remark 4.1. We have considered the situation that a\l\ 7^ 02^2 and a\r\ 7^ 02^2. In case that 
a\l\ = ctili or a.\T\ = a2?"2) then Bl or Br are on the slow manifold Sq and hence there is no 
boundary layer at x = or x = 1. 

4.4 A Special Case 

We conclude the paper by examining a special case. Consider the one-dimensional limit PNP 
system (jl.ip with the special boundary conditions 

0(0) = O , 0(1) = 0, aici(0) = aici(l) = a 2 c 2 (0) = a 2 c 2 (l) = > 0. (4.10) 

One sees that (c±(x) , c^ix) , cj) (x)) = {k/a%, k/d2, (1 — x)4>q) is a steady-state solution. Motivated 
by many works on Lyapunov functions for systems of PNP-types (see, e.g., [3]), we set 
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It turns out that L(t) is a Lyapunov functional. In fact, using the equation and integration by 
parts, 



= ~ M T^KPxCi) H 2~^{o x C2) 

JO V c l c 2 / 



— A / /i(aiCi — «2C2) 2 — A/c / /i(qiCi — a2C2)(ln(a!Ci) — ln(a 2 C2)) < 0. 
io Vo 

Also, L'(i) = if and only if d x c± = d x C2 = ct\C\ — 0202 = 0, which in turn imply that c\ = c®, 
C2 = C2 and (f> = cf)°. Due to the invariant principle (Proposition 12.21 for one-dimensional case), one 
can check that L(t) is equivalent to 



■/ 



(x)) 

if Ci(x, 0) > for x G [0,1], and hence, c, — > c° in L 2 (0, 1) exponentially as i — > cxd. This shows 
a significant difference in asymptotic behavior between the total non-flux boundary conditions 
(see [3]) and the boundary conditions (|3,4p considered in this work for the PNP systems. 
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